clear; close('all'); clc;
addpath '..\common';

fs = 640e6;

iqpath = '.\data\BW200M.wv';
[in_data, in_info] = rm.rfutils.load_rswv(iqpath);
in_data = rm.rfutils.resample_iq(in_data, 600e6, fs); % 原输入波形采样率是 600MHz，重采样到 640MHz
%rm.rfsa.spectrum(in_data, fs, true);

%FS = 625 MHz, Burst 0 startIdx 125986, stopIdx 176506
iqpath = '.\data\RF_Board_2400.iq\File_2025-10-23173051.complex.1ch.float32';
[out_data, out_info] = rm.rfutils.load_iq(iqpath, 'float32', 0, in_info.SAMPLES); % 输出的 IQ 数据时间比较长，只取一个周期的
rm.rfsa.spectrum(out_data, fs, true);

% 计算频率响应
%h_freq = rm.rfsa.ap_response(in_data, out_data);
% 绘图
%rm.rfsa.plot_responses(h_freq, fs);    

% 计算频率响应，并补偿输出幅度和相位
[h_freq, out_compensated]  = rm.hcomp.compensate_rf_response(in_data, out_data, 'smoothing_factor', 0);
% 绘图
rm.rfsa.spectrum(out_compensated, fs, true);
rm.rfsa.plot_responses(h_freq, fs);    
